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<D , Abstract 



Based on an effective Lagrangian which combines chiral SU(3) dynamics 
with vector meson dominance, we have developed a model for the s-wave 
vector meson-nucleon scattering amplitudes. We use this as an input for 
the low energy part of the current-current correlation function in nuclear 
matter. Its spectrum enters directly in the "left hand side" of QCD sum 
rules. For the isovector channel we find a significant enhancement of the 
f ■ >, ' in-medium spectral density below the p resonance while the mass of p 

t*"^ ■ meson itself decreases only slightly. The situation is different in the 

isoscalar channel. Here the mass and therefore the peak position of the 
uj meson decreases strongly while the width increases less drastically 
than in the p meson case. For the meson we find almost no mass 
Qh. shift, the width of the peak broadens due to the imaginary part of the 

(— ! \ scattering length. We find a remarkable degree of consistency with the 

operator product expansion of QCD sum rules in all three channels. 
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H ' 1 Introduction 



At present there is a lively discussion about the behaviour of vector mesons in 
dense and hot hadronic matter. According to Brown and Rho [0 the vector 
meson masses should drop like the pion decay constant in baryonic matter 
(BR-scaling). Other scaling laws || using a bag model find a mass reduction 
of two third times that of the nucleon mass. Several analyses of the p and 
uj meson masses in matter using QCD sum rules 0, |4]] seem to confirm BR- 
scaling. On the other hand, model calculations of the density dependent two 
pion self-energy of the rho meson in nuclear matter || [| suggest only marginal 
changes of the in-medium rho meson mass, but a strongly increased decay 
width instead. Additional BR scaling seems to be needed in order to match 

iWork supported in part by GSI and BMBF 



the hadronic models with the QCD sum rule analysis H. Phenomenologically, 
a dropping p meson mass seems to help ||, (| understanding the enhanced 
dilepton yields seen at masses below the rho meson resonance in the CERES 
and HELIOS experiments at CERN (see the updated review by A. Drees at 
this workshop). 

The QCD sum rule approach, at any practical level so far, has used only a 
caricature of the true isovector spectrum, namely a parameterization in terms 
of a delta function at the meson pole accompanied by a theta function type 
continuum at higher masses. While such a parameterization is valid in vac- 
uum (as we shall confirm), this turns out not to be the case for the in-medium 
spectrum. We will show that using a hadronic model to calculate the vec- 
tor meson-nucleon scattering lengths leads to correlation functions which are 
consistent with the QCD sum rule analysis without need for additional BR 
scaling. 

In section 2 we review properties of the current correlation function in 
vacuum and give a brief account of the QCD sum rule approach. The current- 
current correlation functions in baryonic matter will be developed in section 
3. The comparison with QCD sum rules at finite baryon density is made in 
section 4. 

2 Vector mesons and QCD sum rules in the vacuum 

Before we start looking into the in-medium properties of neutral vector mesons 
it is useful to give a brief reminder of their vacuum properties. Since the vector 
mesons are not stable they can only be seen as resonances in current- current 
(CC) correlation function, 

iWg) =tJd 4 xe^ x (0\TUx)MO)\0) (i) 

where T denotes the time-ordered product and j\i is the electromagnetic cur- 
rent. It can be decomposed as 
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Current conservation implies a transverse tensor structure 
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n(g 2 ), 
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with the scalar CC correlation function 
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U(q 2 ) = -g^IL^q). 



(7) 
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Figure 1: Calculated spectra of current-current correlation functions. The solid 
lines show the vacuum spectra in the p, u and <p channels normalized such that 
they can be compared directly with the corresponding e + e~ — > hadrons data. 
The dashed lines show the spectral functions in nuclear matter at densities 
Po/2 and p — 0.17 fm~ 3 , as discussed in section 3. 



The imaginary part of the correlation function is proportional to the cross 
section for e + e~ — > hadrons: 

„. . cr I=1 (e + e~ — > hadrons) 127r „. . 

R s = — h ^~v^ = Imn s 8 

where a(e + e~ — > p + p~) = AiTa 2 /3s with a = e 2 /47r = 1/137. The vector 
mesons can be distinguished by looking at different hadronic channels with 
corresponding flavour (isospin) quantum numbers. For example, G-parity de- 
mands that the isovector current describing the p meson can only decay into 
even numbers of pions (see data in fig. la). Similarly the annihilation of e + e~ 
into odd numbers of pions determines the isoscalar current describing the uj 
meson. This can be seen in figure lb. Isospin violating processes are small but 
visible as puj mixing corrections to the pion formfactor. The situation for the 
4> meson is more involved. The OZI rule forbids decays into pions, however a 
strong violation of about five percent is seen in the three-pion channel. (Some 
experimental data near the </> resonance have been left out in fig. lb to empha- 
size the cu-contribution given by the solid line). The <fi meson decays mainly 
into OZI allowed channels such as K + K~ or K^Kg (see fig. lc). However 
these channels also exist for the p and uj mesons. The measured cross sec- 
tion of the annihilation into kaons includes an interference between all three 
vector mesons. (For more details see ref fill). Nevertheless the <j) meson still 
dominates the data in these channels. 

Vector meson dominance (VMD) leads to the solid lines in figs. la, b and c. 
In comparison with the experimental data for e + e~ annihilation into hadronic 
channels, we clearly see that the low energy region of the CC correlation func- 



tion is very well described by VMD. The extended VMD model [IT] gives: 
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Imn(g 2 ) = E^V^|^ 2 )r, ( 9 ) 
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where 



Ftf)= (l-°v)g 2 -4 . (10 ) 

q 2 - m v -n^ ac (g 2 ) 

Here my are the bare vector meson masses, n^ c the vacuum self energies and 
gv the strong SU(3) couplings of the vector mesons. The parameter ay is free, 
but usually close to one. It is identical to unity in case of "complete" VMD in 
which all of the hadronic electromagnetic interaction is transmitted through 
vector mesons. For the rho meson channel F(q 2 ) is simply given by the pion 
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where g (g p ) is the coupling of the p meson to the pion (photon) and U p is the 
p meson self energy which is dominated by the two-pion loop [|ll|] . 

In the high energy region the measured correlation function approaches the 
asymptotic plateau predicted by QCD: 
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for large s, where d p = 3/2, d w = 1/6 and d^ = 1/3. 

The basic idea of QCD sum rules is to compare two ways of determining 
Il(g 2 ) in the region of large spacelike Q 2 = —q 2 ^> 1 GeV. One way is to use a 
twice subtracted dispersion relation for each one of the channels V = p, u, <p: 



U(q 2 ) =U(0) + cq 2 + — f ds 
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where the vanishing photon mass in the vacuum requires n(0) = 0. The 
other one is to calculate the correlation function using the operator product 
expansion (OPE): 
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Here the coefficients 01,2,3 incorporate the non-perturbative parts coming from 
the condensates such as the gluon condensate in c<i or the four-quark con- 
densate in C3. The commonly used values 
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channel are c? = d% = 0.04 GeV and C3 = C3 



for the isoscalar and isovector 
-0.07 GeV 6 . The coefficient a 



is proportional to the squared quark mass; due to the small mass of the up and 
down quark we can safely neglect those contributions and set c p { w = 0. The 
situation is different in case of the meson. Because of the large strange quark 
mass, c 2 changes and c\ is no longer negligible. Here we take cf = —0.07 GeV 2 , 
4 = 0.17 GeV 4 and cf = -0.07 GeV 6 PJ. 



A Borel transform is used in order to improve the convergence of the OPE 
series. Comparing eq. ([13|) and (14) after Borel transformation we end up with 
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where the Borel mass parameter should be chosen in the range -M>1 GeV in 
order to ensure convergence of the OPE side of eq.(15). The "left side" comes 



from the dispersion relation of eq.(^) and R represents the ratio (§), but 
now specified for each individual flavour channel with V = p, u>, <fi. The "right 
side" is determined by the OPE. In fig. 3a, b and c we show the comparison 
between the "left side" (solid line) and the "right side" (dashed line) for the 
various (vacuum) channels. The consistency between "left" and "right" sides 
is evidently quite satisfactory. Only for a Borel mass below 0.8 GeV higher 
order terms in the OPE become important and the convergence fails. Keeping 
this success of QCD sum rules in mind we now want to explore whether this 
still holds in the medium. 

3 Current correlation functions in baryonic matter 

The CC-correlation function in medium at temperature T = is defined as 

+oo 

U^(u,q;p) = i f dt [d 3 xe iu)t - i<! -* (16) 

— oo 

* (p|rj M (t,x)^(o)|p), 

where we have replaced the vacuum by \p), the ground state of infinite, isotropic 
and isospin symmetric nuclear matter with density p. We assume matter as 
a whole to be at rest. This specifies the Lorentz frame that we will use in 
the following. For q = 0, the case where the vector meson is at rest, only the 
transverse tensor structure survives and IT 00 = n - 7 = IP = 0. We write 

Uij(uj, q = 0;p) = -8ijn(uj,q= 0; p). (17) 

As a first step we apply the low density theorem to both the correlation func- 
tion and the self energy of the vector mesons. We have 

U{u,q= 0; p) = n vac (w 2 ) - pT{u) + ... , (18) 



and 

IV(cu, q = 0; p) = U^{u 2 ) - pT VN {u) + ... , (19) 

where T(u) is the Compton amplitude and TyN is the vector meson-nucleon 
scattering amplitude, both taken at q — 0. On the other hand using eq.(|19"D 
we can extend eq.(P) to the in-medium correlation function in terms of VMD 
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Comparing the terms linear in density of eqs.(18) and (^0|) we find that the 



Compton amplitude translated into VMD language becomes: 

T(u) = g 2 v F(u 2 ) T VN (u) F(u 2 ). (21) 

In order to determine the in-medium correlation function we are left with the 
task to calculate the vector meson-nucleon scattering amplitudes. We use an 
effective Lagrangian which combines chiral SU(3) dynamics with VMD. This 
Lagrangian has been applied successfully in the vacuum Jll| and has been 
extended to incorporate meson-baryon interactions. For the baryons we include 
nucleons, hyperons and A's. The most important processes contributing to 
the scattering amplitudes are shown in fig. 2 a and b for the p, uj and <fi meson 
respectively. For the p ((f)) meson we draw the diagrams which survive in the 
limit of large baryon mass (fig. 2a-d). The last one (fig. 2d) only contributes to 
the real part of the pN (<f)N) scattering amplitude. While for the p-meson the 
it N and it A loops govern the scattering amplitude, KY* and KK loops dominate 
for the meson. For the uoN scattering amplitude we only show the dominant 
contributions to the imaginary part (figs. 3e, f). We evaluate the imaginary 
parts of the amplitudes by cutting the diagrams in all possible ways. The 
real parts are then determined by using a once subtracted dispersion relation, 
with the subtraction constant fixed by the Thomson limit. Evaluating those 
diagrams and using eq.(^) we plot the spectra of the correlation functions 
for various densities as shown by the dashed curves in fig. 2a, b and c. We 
observe important differences between the various channels. The p meson 
mass decreases only slightly while the width increases very strongly. This 
causes the peak to shift downwards and broaden. We also see strong threshold 
contributions starting at the pion mass. On the other hand the mass of the 
uj meson decreases significantly. Its width also increases but not as strongly 
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Figure 2: Dominant diagrams for the Compton amplitude T(uj) and the vector 
meson-nucleon amplitudes T VN {ui). In order to get from T to T VN , replace 
the (wavy) photon line by the relevant vector mesons. Diagrams (e) and (f) 
operate, in particular, in the u meson channel. 

as for the p meson. For the <fi meson there is almost no change in the peak 
position while the width increases. 

4 Comparison with QCD sum rules 

We now proceed to test the consistency of these in-medium correlation func- 
tions with the OPE of QCD sum rules. Performing an OPE of the density 
dependent correlation functions || the coefficients of eq.(^) become density 
dependent. We use the values d^ = 0.04 + 0.018(p/p ) GeV 4 and d^ u = 
—0.07 + 0.036(p/p )GeV 6 proposed by Hatsuda et al. [Q for the p and ui me- 
son. They neglected higher twist condensates which are hardly known but 
might be important. For the four-quark condensates we assume as in ref.Q 
that ground state saturation holds to the same extent as in the vacuum. For 
the meson the coefficient c\ remains unchanged and we take the values 
4 = 0.17 + 0.01(p/p ) GeV 4 and cf = -0.07 + 0.006(p/p )GeV 6 as suggested 
by Asakawa and Ko fl2|]. In figure 3 we have plotted the Borel transformed 



OPE at normal nuclear matter (long dashed lines) and compared with the "left 
side" (dashed dotted line) as given by the dispersion relation representation of 
our calculated in-medium correlation function. 

The consistency at normal nuclear matter is not as excellent as in the 
vacuum but still very impressive. We therefore conclude that QCD sum rules 
and our hadronic model of the in-medium current-current correlation function 
are mutually compatible. We also point out that assuming a simple pole ansatz 
for the spectrum at finite densities can lead to erroneous interpretation of the 
in-medium masses. This is very obvious in the isovector channel. While our 
model gives only a small change of the rho mass a simple pole ansatz leads to 
a reduction of m p by more than 10 percent || |J at p = p . 
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Figure 3: Comparison of the "left" and "right" side of the QCD sum rules 
(15) in the vacuum (solid and dashed line) and at normal nuclear density (dot 
dashed and long dashed line) as a function of the Borel mass parameter M. 
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